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SUMMARY  AND  CONCLUSIONS 


This  paper  is  concerned  mainly  with  the  "lift/drag"  ratio  (L/D)  which  could  be 
achieved  with  semi -submerged  ships  (S^) .  Given  that  the  motion  in  a  seaway  can 
be  very  small,  then  L/D  is  by  far  the  most  important  factor  in  any  comparison 
with  other  forms  of  vehicle.  That  motions  can  be  less  than  for  conventional 
ships  has  been  clearly  shown  by  a  number  of  eminent  workers  over  the  last 
thirty  years.  We  suggest  here,  in  fact,  that  motions  could  be  essentially 
zero,  if  one  is  prepared  to  accept  a  form  which  has  negligible  hydrostatic 
pitch  and  roll  stiffness  when  underway,  relying  on  dynamic  pressure  forces  on 
stabilizer  foils  for  stability.  The  struts  which  support  the  above-water  por¬ 
tion  of  the  ship  can  then  be  sized  by  structural  considerations  only.  Of 
course,  this  implies  that  the  upper  hull  will  descend  to  the  water  surface  for 
low  speed  operation,  as  indicated  in  Figure  1.  In  many  ways,  operation  of  the 
Figure  1  configuration  would  be  similar  to  operation  of  a  hydrofoil,  and  there 
are  some  interesting  and  potentially  rewarding  trade-offs  to  be  made  between 
buoyant  and  dynamic  lift.  Note  also  that  the  ability  to  employ  a  multiplicity 
of  struts  means  that  the  upper  size  limitations  of  a  conventional  hydrofoil 
are  evaded. 

The  SWATH  configuration  in  Figure  1  has  / 2  times  more  wetted  area  than  the 
minimum  wetted  area  possible,  so  we  show  an  alternative  configuration  in  Figure 
2.  This  was  first  suggested  by  R.W.  Priest  in  the  fifties,*  although  with 
struts  large  enough  to  provide  hydrostatic  stability. 

The  hydrodynamic  efficiency  (L/D)  of  such  a  single  body  is  shown  in  Figures 
3-5,  based  on  calculations  described  later  in  this  paper.  If  we  select 
60  knots  as  the  speed  of  interest,  we  see  that  L/D  increases  markedly  with 
size,  but  also  with  depth  of  immersion  in  the  larger  sizes.  For  A  =  2000  tons, 

L/D  a  10  with  the  body  top  at  the  surface 

a  20  with  the  body  top  48  feet  below  the  surface. 

For  A  =  20,000  tons,  there’s  not  much  difference  when  the  draft  is  shallow, 
but  a  total  draft  of  161  feet  gives  L/D  =  39.  L/D  =  50  is  theoretically 
attainable,  if  the  body  is  deep  enough.  These  figures  are  very  attractive  by 
comparison  with  other  types  of  advanced  vehicle. 

So  far,  we  have  only  considered  conventional  underwater  bodies.  Suppose  now 
we  halve  the  skin  friction  coefficient  in  some  way;  either  by  polymer  injec¬ 
tion**  or  by  appropriately  shaping  the  body.  The  ultra- low  drag  underwater 
bodies+are  an  example  of  such  drag  reduction  by  shaping,  by  virtue  of  extensive 

♦Reported  by  Boericke*. 
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**  See  Van  Mater  for  example. 

t  Payne^  describes  the  antecedants  of  this  technology. 
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Figure  1.  A  Small  Waterplane  Twin  Hull  (SltfATH)  Configuration  Which  Becomes 
"Foil -Borne”  for  High  Speed  Operation. 
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Figure  3.  Performance  of  a  200  Ton  Submerged  Body. 
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Figure  5.  Performance  of  a  20,000  Ton  Submerged  Body. 
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laminar  flow  in  their  boundary  layers.  Perhaps  a  more  practical  possibility 
for  S3  is  to  design  for  a  "tired,"  low  friction  turbulent  boundary  layer  over 
most  of  the  body.  Stratford^  has  demonstrated  such  a  flow  experimentally,  and 
while  no  one  has  yet  investigated  the  technology  exhaustively  for  external 
flows,  it  would  appear  to  offer  some  promise. 


A  third  alternative  is  drag  reduction  by  injected  air  lubrication. 

If  such  a  halving  of  the  skin  friction  coefficient  is  achieved,  then  the  L/D 
ratios  for  the  three  different  displacements  become  as  shown  in  Figures  6-8. 
Even  the  200  ton  size  is  now  competitive  with  other  vehicles  at  60  knots 
(L/D  -  17)  while  values  as  high  as  L/D  =  80  are  feasible  in  the  20,000  ton  size. 


It  should  be  emphasized  that  all  of  these  values  will  be  degraded  by 

•  Strut  drag. 

•  Foil  drag. 

•  Resistance  of  any  other  appendages. 

By  careful  design,  however,  we  can  minimize  the  penalties  involved.  But  since 
strut  size  dominates  strut  drag,  and  strut  size  depends  principally  on  struc¬ 
tural  loads,  we  have  not  attempted  to  include  estimates  for  these  parasitic 
items  in  this  paper. 

3 

We  conclude  that  in  medium  and  large  sizes,  S  offers  the  possibility  of  very 
high  L/D  ratios  in  the  intermediate  speed  range  around  60  knots,  and  large 
S3  may  be  competitive  at  100  knots.  The  poor  showing  to  date  seems  to  be 
principally  due  to  designs  which  have  the  submerged  hull  too  close  to  the 
surface,  struts  which  are  very  thick  and/or  have  large  wetted  area,  and 
insufficient  development  work  on  reducing  interference  drag. 
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Figure  8.  Performance  of  a  20,000  Ton  Submerged  Body  if  the  Skin  Friction 
Coefficient  is  Halved. 
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A  VERY  BRIEF  HISTORY  OF  SEMI-SUBMERGED  SHIPS 


There  does  not  seem  to  be  any  totally  satisfactory  name  for  what  we  are  dis¬ 
cussing.  The  abbreviation  S^  does  at  least  have  the  merit  of  extreme  brevity, 
and  so  we  employ  it  here.  But  we  are  really  talking  about  a  ship  which  has  low 
self-making  wave  drag,  and  a  reduced  response  to  the  seaway.  A  "small  water- 
plane  area"  ship  will  generally  have  both  these  virtues,  but  the  description 
is  very  imprecise.  Also  there  is  at  least  one  other  way  of  achieving  the  same 
result,  as  we  shall  see. 

3 

The  first  patent  for  S  was  awarded  to  Reuben  N.  Perley  (1922)  for  a  submerged 
monohull  with  surface-piercing  struts  forward  and  aft.  A  similar  ship  with 
a  single  strut  (the  "Shark  form")  was  patented  by  Rudolf  Engleman  in  1937.  A 
modern  SWATH  configuration  was  patented  (in  England)  by  Frederick  G.  Creed  in 
1944.  Thus  there  is  nothing  particularly  "new"  about  the  basic  idea  or  concept. 
The  problem  has  been  in  assessing  its  virtues  and  vices. 

It  would  seem  that  the  earliest  studies  in  this  country  were  carried  out  by  a 
group  in  BuShips  (Code  420)  under  Owen  Oakley,  in  the  early  fifties.  Boericke* 
summarizes  much  of  the  work  done  by  this  unusually  creative  team,  which  included 
Robert  W.  Priest,  David  Winter,  James  L.  Mills  and  others.  Unfortunately,  they 
were  apparently  unable  to  arouse  any  interest  in  supporting  the  broad-based 
R^D  effort  which  would  have  been  required  to  select  the  best  among  their  various 
proposed  configurations  and  then  to  develop  it.  Figure  9  depicts  some  of  the 
configurations  studied  at  this  time. 

From  a  different  vantage  point,  Lewis’  studies5  of  seakeeping  in  the  mid-fif¬ 
ties  led  to  the  first  semi -submerged  ship  (S^)  described  as  such,  and  a  form 
which  he  patented  in  1959. ^  The  lines  are  given  in  Figure  10.  This  was  the 
first  time,  apparently,  that  a  form  had  been  consciously  designed  for  super¬ 
critical  operation  in  head  seas;  although  supercritical  operation  was  known 
to  be  possible  with  conventional  forms  under  certain  rather  special  conditions. 
Insofar  as  one  can  tie  it  down  to  one  individual,  Lewis  was  clearly  the  inven¬ 
tor  of  the  supercritical  ship. 

The  work  of  this  intellectually  fruitful  decade  was  summarized  in  a  number  of 
important  papers  which  appeared  in  the  period  1959-1962.  Mandel’s7  study,  and 
that  of  Lewis  and  Breslin°  are  the  best  and  the  most  complete.  But  Mandel’s 
statement  that  "none  of  the  new  ship  types  ...  is  likely  to  supplant  more 
traditional  ships  and  aircraft  .  .  ."  may  have  discouraged  some  further  work, 
because  it  was  so  easily  taken  out  of  context. 

At  this  point,  one  might  say  that  the  most  severe  seakeeping  problem  -  high 
speed  in  head  seas  -  had  been  solved  by  Lewis’  invention  of  the  supercritical 
ship  (SCS)  and  that  this  concept  could  be  applied  to  some  of  the  configurations 
already  proposed.  But  the  monohull  SCS  was  very  tender  in  roll,  had  insuffi¬ 
cient  deck  space  for  many  purposes,  and  experienced  rather  extreme  motions 
at  resonance. 
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Figure  9.  Configurations  Studied  by  the  Code  420  (Preliminary  Design)  Team 
in  the  Early  Fifties.  (From  Bocricke1) 


WLP  *  24  (4.30)  DECK 


Figure  10.  Lines  of  Lewis'  Semi -Submerged  Ship. 


In  the  period  1964-66,  Payne  built  two  small,  manned,  supercritical  catamarans 
(FICAT  I  and  II)  which  avoided  the  first  two  of  these  three  problems.  The 
first  of  these  craft  is  shown  in  Figure  11.  A  reduction  in  wave  drag  was 
achieved  partly  by  the  reduction  in  individual  hull  beam  which  is  permitted 
by  the  catamaran  configuration,  and  partly  by  cancellation  interference  of 
the  wave  trains  from  each  hull.  In  other  words,  interference  was  substituted 
for  submergence  of  the  main  buoyancy  volume.  This  avoided  some  of  the  practical 
problems  associated  with  struts,  excessive  draft,  etc.,  and  resulted  in  a 
more  conventional  ship. 

The  two  FICAT' s  were  operated  "at  sea"  in  the  Chesapeake  area  for  some  hundreds 
of  hours,  during  the  years  1965-66,  under  widely  varying  conditions.  Model 
tests  (Figure  12)  were  also  conducted,  and  it  seemed  clear  that  at  least  a  50% 
wave  drag  cancellation  was  being  achieved  for  all  Froude  numbers  above  about 
0.5,  with  some  evidence  of  total  cancellation  at  F  =  0.7.®  The  FICAT's  were 
the  first  SCS  configurations  to  be  analyzed  theoretically1®* 11  > So  far  as 
resistance  was  concerned.  Band  found  generally  good  agreement  between  Michell's 
wave  drag  integral  and  experiment,  including  wave  drag  calculated  from  photo¬ 
graphs  of  the  local  surface  elevation. 

Despite  the  significant  progress  made,  the  Payne  team  was  never  able  to  find 
support  for  its  research,  which  was  therefore  discontinued  in  early  1969,  so 
that  attention  could  be  concentrated  on  supercritical  planing  hulls.  A 
definitive  assessment  of  FICAT  vis-a-vis  other  configurations  was  never  made, 
and  so  this  must  remain  a  question  mark. 

Coincident  with  Payne's  departure  from  the  field,  Leopold13,14  re-introduced 
the  Creed  configuration  and  built  the  first  Litton  TRISEC15  man-carrying  model. 
This  was  the  first  true  S5  actually  reduced  to  practice,  and  its  seakeeping 
ability  was  very  impressive.  The  problem  of  excessive  motions  near  resonance 
was  solved  by  using  inclined  active  foils  which  acted  as  both  rudders  and 
pitch  angle  control. 

Litton  was  also  unsuccessful  in  obtaining  funding  to  pursue  the  research,  but 
the  basic  concept  has  since  been  an  ongoing  project  in  the  Navy,  principally 
at  NSRDC  and  NAVSEC. 

1  17 

Shortly  after  Leopold,  Lang  *  patented  some  very  innovative  improvements, 
and  was  later  able  to  build  a  190  ton  "manned  model."  Again,  because  of  lack 
of  funds,  this  design  was  not  "optimized,"  and  has  a  much  higher  resistance 
than  the  minimum  possible.  Lang  adopted  and  improved  the  active  foil  stabili¬ 
zation  concept,  and  linked  it  to  a  simple  autopilot  to  achieve  minimum  motions 
at  all  speeds. 

In  this  present  paper,  we  have  suggested  that,  since  active  stabilization  is 
essential  while  underway,  one  might  as  well  dispense  with  hydrostatic  stability 
and  thus  avoid  the  rather  severe  resistance  penalties  which  it  imposes.  We  also 
suggest  that  the  hulls  should  be  much  deeper  in  the  water  if  high  efficiency 
is  to  be  achieved.  Then,  as  a  final  improvement,  the  operational  limitations 
imposed  by  this  deep  draft  can  be  ameliorated  by  hinging  the  strut  assembly  in 
such  a  way  that  the  above  and  underwater  hulls  can  be  brought  together  for 
low  speed,  low  draft  operation.  Although  there  may  be  some  size  limitations 
to  this  last  suggestion,  the  basic  technology  is  not  far  removed  from  that  used 
in  the  new  variable  wing  sweep  bombers. 
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THE  PERFORMANCE  CALCULATIONS 


Oddly  enough,  the  basic  theoretical  tools  for  an  overview  assessment  were  devel¬ 
oped  by  Havelock  prior  to  the  need  for  them.  The  four  relevant  papers  are: 

18 

’’The  Wave  Resistance  of  a  Spheroid" 

19 

"The  Wave  Resistance  of  an  Ellipsoid" 

20 

"The  Moment  on  a  Submerged  Body  Moving  Horizontally" 

21 

"The  Forces  on  a  Submerged  Body  Moving  Under  Waves" 

For  the  present  paper,  we  are  interested  in  wave  resistance  only,  the  equations 
for  which  are  given  in  Appendices  I  and  II.  Since  the  step  from  Havelock's 
equations  to  numerical  results  is  not  entirely  trivial,  we  first  checked  our 
results  with  the  only  known  previous  solution,  the  "slender  body"  numerical 
evaluation  by  Wigley.^2  As  Figure  14  shows,  the  agreement  is  good. 

As  recounted  in  Appendix  I,  we  then  compared  the  theory  with  tank  test  measure¬ 
ments  of  residual  resistance;  again  with  satisfactory  agreement  for  our  present 
purposes. 

Some  general  trends  were  then  studied.  Figure  15  gives  numerical  results  for  a 
body  of  revolution,  and  Figure  16  shows  the  effect  of  varying  the  cross-sectional 
shape  for  a  fixed  submergence  of  the  centerline.  The  corresponding  total  resis¬ 
tance  ratio  is  given  in  Figure  17,  and  the  inverse,  L/D  in  Figure  18.  Changing 
the  cross-sectional  shape  clearly  has  a  very  small  effect  compared  with  a 
change  in  depth,  so  this  effect  was  ignored  in  subsequent  calculations. 

The  final  calculations  are  summarized  in  Figures  3-8,  which  have  already  been 
discussed  at  the  beginning  of  this  note. 

Some  limited  work  on  strut  resistance  was  done  during  the  course  of  this  work, 
and  is  summarized  in  Appendix  III.  This  is  considered  to  be  incomplete. 
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Figure  14.  Wigley's  "Slender  Body"  Numerical  Approximation  in  Comparison 
with  the  Exact  Solution  of  Havelock's  Equation  for  the  Wave 
Drag  of  a  Prolate  Ellipsoid. 
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Figure  15.  The  Effect  of  Slenderness  and  Submergence  on  the  Wave  Drag 
of  a  Prolate  Ellipsoid. 
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Figure  16.  The  Effect  of  Varying  the  Cross-Secti 
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APPENDIX  I 


CALCULATION  OF  SMALL  WATER  PLANE  HULL  PERFORMANCE 
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In  considering  the  feasibility  of  a  new  ship  concept,  an  early  question  is 
"what  is  its  transport  efficiency?"  We  typically  evaluate  the  calm  water 
"lift/drag  ratio"  (the  inverse  of  the  more  conventional  R/D)  in  order  to 
answer  this,  and  it's  clear  that  approximate  figures  for  "L/D"  are  quite 
adequate  for  establishing  feasibility  so  long  as  they  are  realistic.  In  the 
case  of  S3,  it  is  rather  simple  to  compute  resistance,  and  it's  also  clear 
that,  due  to  small  motions  in  a  seaway,  performance  will  not  degrade  much 
with  increasing  sea  state;  except  possibly  near  i*esonance. 


The  significant  resistance  components  are 


Rw  Hull  wave  making  drag 

°BS  | HuU  pressure  draf  \  (considered  together  in  this  analysis) 

Strut  wave  making  and  spray  drag 

D  | Strut  skin  friction  drag 

SS  |  Strut  pressure  drag 

Dj  Strut/hull  interference  drag 

D.  Wind  resistance 

A 

The  purpose  of  this  Appendix  is  to  present  equations  for  the  first  three  resis¬ 
tance  components. 


Hull  Wavemaking  Drag 


This  is  perhaps  the  most  critical  term,  because  it  will  dominate  the  optimization 
of  strut  length.*  On  the  other  hand,  it  would  be  needlessly  complicated  to 
determine  the  wave  drag  of  different  hull  shapes  at  this  stage  in  the  analysis, 
when  we  are  looking  for  overall  trends.  It's  sufficient  to  determine  a 
"standard"  variation  of  wave  drag  with  slenderness  ratio  and  immersion  depth, 
recognizing  that  subsequent  variations  of  hull  volume  distribution  may  enable  us 
to  improve  on  this  result. 

Such  a  "standard"  variation  is  provided  by  Havelock's  analysis  of  a  prolate 
spheroid1’1,  where  he  obtains 

.  128,2gpaWe-2'f2/ e-2'2/f2  [JV1WI2  ■"  « 

o 


*  A  trade-off  which  does  not  seem  to  have  been  addressed  in  the  literature. 
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where 


e  =  Y\  -  (b/a)2  the  eccentricity 

a,b  are  the  major  and  minor  semi-axis  lengths 

f  =  u/t/gh,  the  submergence  Froude  number 
h  =  submergence  of  the  spheroid  centerline 
^2/2(2)  is  the  Bessel  Function  of  the  first  kind,  of  order  3/2 
=  V2/irz  -  cos  zj 


t  is  a  dummy  variable. 


The  displacement  of  a  prolate  ellipsoid  is 


A  = 


(1.2) 


Thus,  from  equations  1.1  and  1.2 

RW  |  / a  \2  3,2(  -2/f2  2  [  e'2t  /f  |"sin  z  l2  w  „ 

T  =  j  967r(b)EA  je  .  cos  zj  dt  (1.3) 

I  )  0 

The  term  in  the  |  j  bracket  is  a  function  only  of  b/a  and  is  plotted  in  Figure 
1.1,  and  tabulated  in  Table  1.1. 

Equation  1.3  is  easy  to  integrate  numerically  as  it  stands,  so  that  there  is 
little  point  in  seeking  further  simplification.*  In  Figure  1.2,  we  have  com¬ 
pared  it  with  some  experimental  measurements  of  net  residual  resistance  from 
Appendix  7  of  Reference  1.2.  Since  the  experimental  data  is  for  "streamline" 
bodies,  the  aft  portions  of  which  are  quite  unlike  an  ellipsoid,  the  agreement 
seems  excellent,  and  quite  sufficient  for  our  present  purposes. 


*  Because  of  the  exponential  term,  it  converges  well  as  t  increases. 
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Table  1.1. 


Functions  of  the  Fineness  Ratio  b/a 


(b/a) 

e 

A 

{f (b/a)  } 

IE-03 

. 9999995 

2 . 500021 IE-07 

1 . 8849846E 

IE-02 

.99995 

2. 5014505E-05 

1 . 8868605E 

5E-02 

.99874922 

6. 3161389E-04 

4.7946157E 

IE-01 

.99498744 

2. 5905254E-03 

.19936531 

.2 

.9797959 

1. 1260011E-02 

.89917777 

.3 

.9539392 

2. 8651639E-02 

2.3880292 

.4 

.91651514 

6.0075642E-02 

5.2374138 

.5 

. 8660254 

.11643375 

10.622599 

.6 

.8 

.22249714 

21.234275 

.7 

.71414284 

.44496157 

44.383883 

.8 

.6 

1.0231108 

106.54675 

.9 

.43588989 

3.5216047 

382.42625 

.95 

.3122499 

10.894751 

1207.5788 

.98 

.19899749 

45.337269 

5086.5397 

.99 

.14106736 

130.40137 

14689.081 

MOOEL  4157 

SUBMERGENCE  DEPTH  RATIO 

- O  h/b=  2.544 

- □h/b  =  2.84 

- A  ^  =3.164 

3. 


THEORY 


^EXPERIMENT 


<  O.OI 


LI 


5  «  7~1  '  (0  2  3  4  8  6  7  «  9  |0.0 

SUBMERGENCE  FROUDE  NUMBER  Fh=u///giT 


Figure  1.2.  Comparison  Between  the  Theoretical  Wave  Drag  of  a  7:1  Prolate 
Ellipsoid  and  Some  Experimental  Measurements. 


WAVE-MAKING  RESISTANCE  RATIO 


Figure  1.2.  Continued. 
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Hull  Friction  and  Pressure  Drag 
I  3 

Hoerner  gives  a  relationship  for  this  which  has  been  widely  accepted.  De¬ 
fining 


Dn„  =  Cn  i  pu2  S  „ 

BS  D  2  wet 

wet 


:fj1  +  l(!)5/2*7(!)s] 


(1.4) 


(1.5) 


where  is  the  usual  flat  plate  skin  friction  coefficient,  which  may  be  con- 
B  veniently  expressed  as 


_ 0.427 _ 

fB  dogl0R2a  -  0.407)2*64 


(1.6) 


R2a  =  — —  (Reynolds  number  based  on  length) 


(1.7) 


-S  2 

v  «  1.25  x  10  (ft  /sec)  for  normal  water 


There  is  no  simple  relationship  between  the  wetted  area  S  r  t  and  b/a,  because 
other  factors  are  involved,  principally  the  prismatic  coefficient. 


Volume 

„  ,2 


(1.8) 


A  typical  variation  of 


(1.9) 


with  b/a  and  C„  as  given  in  Figure  1.3.  Cg  is  obviously  less  at  the  lower 
prismatics;  bur  then,  so  is  the  displacement,  so  this  tends  to  cancel  out. 
For  the  purpose  of  preliminary  performance  calculations,  therefore,  we  might 
as  well  use  the  spheroid  relationship 


S  ^  =  2*a2  (*)(*  + 

wet  \  a  "a  e  / 


(I. 10) 


where 


(b/a) 4 


Since  we  are  generally  working  to  a  given  displacement,  the  hull  volume  V  is 
known,  and  since 

•  •  isr  (S)i/s 
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O  0.65 


0  0.7 


DATA  FOR  SERIES  58  FORMS  FROM  REF.  i.2 


PRISMATIC  COEFFICIENT  Cp  =  VOL /27r  ab2 


JNDER  WITH  SPHERICAL  ENDS 

CP='-T'T> 


"SPHERE 


PROLATE  SPHERIOD 
CP  =  2/3^\ 


•SLENDER  QUADRATIC  SPINDLE 


C  =-§- 
CP  15 


DIAMETER  _  b. 
LENGTH  "  a 


Figure  1.3.  Wetted  Area  Coefficient,  as  a  Function  of  Fineness  Ratio 
and  Prismatic  Coefficient. 
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APPENDIX  II 


WAVE  RESISTANCE  OF  SUBMERGED  ELLIPSOIDS 
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The  problem  of  the  wave  resistance  of  generalized  ellipsoids  has  been  solved  by 
Havelock** *1.  The  purpose  of  this  Appendix  is  to  recast  Ids  equations  into  a 
form  suitable  for  numerical  evaluation. 

Case  of  a  >  b  >  c 


From  Havelock 


2  2  2  2 
„  _  32n  gpaVc 


32lr  g£2  -o— 2  x/2  e'2Koh[Int  (1)  +  Int  (2)] 
(2-a  (a-oy'  ^ 


(II. 1) 


4 

The  volume  of  the  ellipsoid  is  j  irabc.  Thus 


24rrabc _ 

r9  i2/-.2  k2a3/2 
(2-c*o)  (a  -b  ) 


e"2K°h  [Int  (1)  +  Int  (2)] 


(11.2) 


In  Havelock’s  notation  K  =  g/u2.  We  will  employ  the  submergence  Froude  number 
f  =  u//gh,  so  that 


o  r2  . 
f  gh 


2KQh  =  2/f 


(II. 3) 


For  numerical  integration 

f  A 


f i 

"bc/  ?;  2  ;,3 
<  /(a  +a) 


dA  . _  .  2  .  ,  -3/2 

2  2  +  3  a*,cAi 

(b^+A) (c  +A) 


(11.4) 


2  2  2  .  1 1  2 
where  Aj  >>  a  ,b  ,c  .  (This  form  of  the  well-known  jntegral  is  due  to  Band  .) 

In  nondimcnsional  form ,  if  b  =  b/a  c  =  c/a  A  =  A/a 


a  =  be 
o 


/  dA 

J  /(l+A)J(b2+X)(c2*A] 


2  ;  -3/2 

+  -bcAi 


(11.5) 


=  - 2f-~y  wTT-  c-2/f'  [Int  (1)  ♦  Int  (2)] 

&  (2-ao)2(l-b2)V2 


(1 1*6) 
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Defining 


1 


-.2  2  r-2  -2 

2  b  -c  b  -c 


1  *2  Vi2 

a  -b 


-2 

l-bz 


(II. 7) 


Havelock  gives 


Int(l)  = 


1/a 

■/ 


‘[J3/2  j^2 


1  ,-2  J2 


2.2,,l/2l2  e'2t2/f2  dt 


(1  -  ai2t2)3/2 


where  the  Bessel  Function  J.,„  is  defined  in  Appendix  I.  Thus  the  integral  can 
be  written  as  ' 

l/“l 

■»'2t  l(si"  dt  (II. 8) 

<1  (l-c.12t2)3/2 


Intel)  =  i 


1/u 

■1/ 


where 


(l-b2)(l+t2)(l-ai2t2) 

7?” 


(II. 9) 


(h  =  h/a) 

For  the  second  integral,  Havelock  gives 

,2  -2t2/f2 

(P)l  e 


Int  (2)  = 

1/a 


f  [!j/2  (p 

/  (a2t 


1 


,  ,.2  ..3/2 

(aj  t  -1) 


dt 


(II. 10) 


where 


■  V 


(l-b2)(l+t2)(a12t2-l) 

7? 


(ii. id 


and  I3/2  is  the  modified  Bessel  Function  of  the  first  kind,  of  order  3/2 
Since 

Ij^2(z)  =  ^2/ttz  sinh  z,  I  |^2(z)  =  /2/irz  cosh  z 


and 


2n 


I  ,  =  1  .  -  —  I 

n+1  n-1  z  n 


I  -I  .  I  1 
*3/2  "  -1/2  z  1/2 


(a  standard  form) 
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I3/2(-P')  =  [cosh  (p) 


sinh(p) 


] 


(II. 12) 


Int  (2)  = 


co 

1/a, 


,-2t2/f2  [cosh(p)  -  S1— p-^-  ]2  dt 


(a,2t2-i)s'2 


(11.13) 


Equation  II. 6  can  now  be  evaluated,  using  equations  II. 5,  II. 8  and  11.13. 


Case  of  a  >  c  >  b 


From  Havelock,  after  some  manipulation 


where 


_W 

A 


Int (3)  = 


24vbc  e~2/f 
^2-ao)2(l-b2)3/2 


Int (3) 


tJ3/2(r)1  6 


/  O  4  “22t2) 


2  -2t2/f2 


dt 


where 


and 


2  f  e-2t 

~  -nj  r 


-,J2.,C2  / sin  r 

_2t  /f  ( — - - cos  r | 


(1  .  a22t2)3'2 


dt 


1  -  b2 


r  _./(l-b2)(ln2)(haV 

V  f4h2 


Thus  IT. 14  can  be  evaluated  using  11.15  and  11.16. 


(11.14) 


(11.15) 


(11.16) 
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Wave  Drag  of  a  Sphere 


From  Havelock,  after  some  manipulation 

00 

R  _  3  -2/f2  f  n  .2.3/2  -2t2/f2  , 

1  ~  ?  /  (1  5  6  dt  (II. 

o 

where  f  =  u/»^gh,  the  submergence  Froude  number. 

This  solution  is  a  useful  check  case  for  the  general  equations  when  programmed 
for  numerical  solution. 


APPENDIX  II  REFERENCES 


11. 1  Havelock,  T.H. 

11. 2  Band,  E.G.U.,  and 
Payne,  P.R., 


"The  Wave  Resistance  of  an  Ellipsoid."  Proc. 
of  the  Royal  Society,  A.,  Vol.  132  (1931). 

"The  Pressure  Distribution  on  the  Surface  of 
an  Ellipsoid  in  Inviscid  Flow."  Payne  Inc. 
Working  Paper  No.  101-15  (July  1972). 


42 


APPENDIX  III 

WEDGE  STRUT  RESISTANCE  ESTIMATE 
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During  this  study,  it  was  felt  that  a  base-ventilated  wedge  section  strut  might 
have  some  attractions.  Accordingly,  the  following  calculation  was  carried  out. 


Surface  Pressure  Drag 

The  equations  for  the  force  per  unit  length,  P/J,,  on  a  wedge  are  given  by 
Korvin-Kroukovskv  and  Chabrow* I I • 1 ;  viz 


where 


P/Jt  =  2pkbv  cos  g  f(y) 


(III.l) 


n  -  1  _  £ 

P  o  n 


f(Y) 


1 

k 


2  2  and  0  as  t^ie  wec*ge  included  angle 

ir/2 


_  f  f/1  +  sin  y  \  _  / 1  +  sin  y\  1 

J  cos  Y  /  \  cos  Y  /  J 

o 

//2 

a* 


cos  y  sin  y  dY 


4  cos  g  /  (1  +  sin  Y)n(cos  Y)1_n  sin  Y  dY 


n  = 


ir  -  2g 


o_ 

ir 


b  is  the  total  strut  width 

Figure  III.l  and  Table  III.l  give  the  variation  of  P  with  0  obtained  from 
numerical  integration  of  these  equations. 

Some  curve  fits  to  these  results  are  of  value. 

For  0  <  0  <  30° 

c  _  _ _ 

P  "  (89.9543  +  0.517358  0°) 

For  0  <  0  <  180° 

Cp  =  1.1107  x  10-2  (0°)  -  6.1332  x  1O'5(0°)2 

+  2.31799  x  1O"7(00)3  -  6.02678  x  1O"1O(00)4 
+  7.85803  x  10"13  (0 ° ) 5 
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Table  III.l 


(deg.)  Cp 


2 

0.021978 

4 

0.043473 

6 

0.064498 

8 

0.085062 

10 

0.105178 

12 

0.124855 

14 

0.144104 

16 

0.162935 

18 

0.181357 

20 

0.199379 

22 

0.217011 

24 

0.234262 

26 

0.251140 

28 

0.267656 

30 

0.283816 

40 

0.359552 

50 

0.427527 

60 

0.488569 

70 

0.543405 

80 

0.592685 

90 

0.636976 

100 

0.676800 

110 

0.712603 

120 

0.744788 

130 

0.773720 

140 

0.799723 

150 

0.823082 

160 

0.844058 

170 

0.862884 

180 

0.879762 
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Surface  Skin  Friction 

Korvin-Kroukovsky  and  Chabrow***'*  give  the  surface  velocity  as 


/  cos  y  \  n 
\ 1  +  sin  y  / 


where  x  is  related  to  the  dummy  variable  y  by 

/*/2 


(III. 2) 


x  =  2kb  cos  g 


f  ■ 


sin  y)n(cos  y)*’n  sin  y  dy 


fill. 3) 


In  principle,  then,  we  can  compute  the  local  skin  friction  force  and  integrate 
it.  For  the  present  study,  because  of  time  limitations,  we  employ  an  approxi¬ 
mation  which  is  simpler,  although  less  accurate.  We  assume  that  the  surface 
pressure  pg  is  a  constant  so  that  the  pressure  drag  P  is  given  by 


Since 


b£  Ps  P» 


1  2  12 

P„  +  T  Pu  =  P_  +  T  Puc 
00  2  s  2  s 


(III .4) 


w 


Poo  -  ps  +  \  p"2 


1  2  P 

"  2  PU  '  b£ 


2  pu  t1  '  cpl  where  Cp  - 

2  pu  b£ 

therefore  total  friction  drag 

Df  «  (2c£)  Cf  i  pu2  [1  -  Cp] 


(III. 5) 


(III. 6) 


where  c  is  the  chord  and  i  the  wetted  length. 
Based  on  frontal  area,  since 


2  tan  (0/2) 


“  ipu2b£ 


tan  (0/2)  (1  '  CP) 


(III. 7) 


Base  Drag 


Figure  III. 2.  Side  View  of  Strut. 


IVe  can  compute  the  base  pressure  Apg  using  lioerner's  approximation,  or  directly 
from  experimental  observations.  From  Hocrner's  approximation  (p.  3-21) 


Then 


C 


f 


2  tan  (9/2)  cos  (0/2) 


0.135 


tan  (6/2)  cos  (8/2) 

C, 


(1 1 1. 8) 


(III. 9) 


As  Figure  III. 3  shows,  this  agrees  well  with  experimental  observations. 
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SE  PRESSURE  COEFFICIENT  C p^/\p 


The  depth  z  of  the  ventilation  is  then  given  by 


i.e. 


APB  =  Pgz 


zg_ 


=  0.0675 


2  tan  (9/2)  cos 


(9/2)  j 


1/3 


(III. 10) 


This  is  plotted  as  an  equivalent  Froude  number  in  Figure  III. 4. 

u  //gz  is  a  kind  of  Froude  number,  therefore,  dependent  only  upon  0  and  C^. 
Typical  values  are: 

For  Cf  =  .004,  6/2  =  5°  and 


u  = 

10 

20 

30 

40 

50 

knots 

z  = 

2.1 

8.4 

18.9 

33.6 

52.5 

ft 

The  total  base  drag  of  the  strut  will  be 

db  =  {\  2Pg)bz  +  b(fc-z)  ApR 

=  \  Pgbz2  +  b(«,-r)  j  pu2 

1  *  2xT 

=  pgbz2  +  j  pgbz(£-z) 

=  i  pgb(z2  +  £z)  (III. 11) 

where 

z  is  given  by  equation  I II. 10. 
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Spray  Drag 

I II .  2 

From  Ogilvie  ‘  we  take  the  spray  surface  elevation  to  be  (see  Figure  I II. 2) 


-  I  j  cos  (flxX/2)  ^  ]  dX 


(III. 12) 


where  K  =  g/u 

The  wetted  area  associated  with  this  is 

r 

AS  =  /  £  dX 


(III. 13) 


In  the  numerical  integration  of  equation  III. 12,  a  singularity  occurs  at  A  =  0. 
This  difficulty  can  be  avoided  by  integrating  to  infinity  from  a  small  value 
Aj,  and  determining  the  residue  as  follows.  Since  A  <<  1,  cos  (OxA/2)  =  1.0, 
except  right  at  the  bow  (x  =  0) .  Also  sin  (/i(A  x)  =  /iCA  x. 

.1 


At 


■  ?/(— X--) 


ex  /  /  HA 

fO 

/—n 

*< 

CM 

/ — X 

w 

t 

(HA)4 

*  J  \  1! 

2! 

3! 

4! 

o 

L  - 

(HAj)2 

-  --  -  -  j 

(HAj)3 

(HAj)4 

*  l  1 

2.2! 

3.3! 

4.4! 

(III. 14) 


Equation  111,12  is  plotted  in  Figure  III. 5  for  some  arbitrary  values  of  II  and  x. 
Total  Strut  Resistance 

If  z  <  £  (or  gz/u2  <  g£/u2  ,  or  u2/gz  >  uVgfc) 

D  =  2 pl'2  b*cp  + 1 pu2  cf  ?inrre7Tf (1  ■  cp 


2  — p  .  2  "f  tarTT07Tf  ^  ‘  lP 

+  |jpgz2b  +  b(£-z)  ~  pu2  (zg/2u2)  |  +  (spray  drag)(III.15 
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WAVE  HEIGHT 


DISTANCE  ALONG  WEDGE  IN  FEET  -  X 


Figure  III. 5.  Spray  Sheet  Surface  Elevation  from  Equation  111.12. 
t,  =  Praft  in  feet 

This  case  is  for  u  =  25  ft/sec,  O  =  7.5° 
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Some  typical  results  are  given  in  Figure  1 1 1.6.  Spray  drag  has  been  omitted 
from  the  calculation  because  of  time  limitations. 
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Figure  III. 6.  Total  Resist 
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